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The textbook starts immediately with the description of bar and beam elements. We will first 
do an even simpler element — the spring element. The reason for this is that we already know 
the stiffness of the spring k. With several springs in a system, the global stiffness matrix 
needs to be assembled. This procedure can easily be explained using spring elements. We 
will then move onto the bar and beam elements where we will calculate the stiffness matrix 
for each respectively. The process of assembling the element stiffness matrices into a global 


stiffness matrix is the same for any type of element. 


B.1 Spring Elements and Assembly 


B.1.1 Single Element 








Figure B.1-1— Single spring element with stiffness k 


Assume that the spring is linear elastic, following the relation F =kA. Let node 2 be 
displaced by a distance u2 and node 1 by a distance uw; with u2 > u1. Consider the equilibrium 


of forces for the spring under these displacements. 
Ji =k, —u,) 
fa = k(u, — u) B.1-1 


Equation B.1-1 can be written in matrix notation, 
k =k || u f 
parm B. 1-2 
-k k jiu, h 
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stiffness matrix x displacement vector = force vector 


Lets analyse Equation B.1-2 by setting u2 = 0, 


SRH > EHSF ps 


and similarly by setting u = 0, 


bek EL pra 


From this it can be seen that the first column in the stiffness matrix gives the force at each 
node in order to have a unit displacement at node 1 and zero displacement at all the other 
nodes. Similarly, the second column in the stiffness matrix gives the nodal forces with a unit 
displacement at node 2 and zero displacement at all the other nodes. This procedure will be 


used later to determine the stiffness matrix for beam elements. 


3 of 27 


Eindige Element Metodes 414 / Finite Element Methods 414 


B.1.2 System of Spring Elements 
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Figure B.1-2 — System of single spring elements 


For element 1, the element nodes are 1 and 2 as in Figure B.1-1. These nodes correspond to 
the global nodes 1 and 2, as indicated in Figure B.1-2. For this element, the equilibrium 


equation is given by 


k = k u elel 
| | l- J B.1-4a 
-k k fu% h 
For element 2, the element nodes 1 and 2 correspond to global nodes 2 and 3 respectively, 
k ve k u ele2 
h se a 
-k, k, ||“ fh 


Equilibrium of nodal forces can be written as 


elel elel ele2 ele2 


F=f, s F=f th s B=f B.1-5 
where F indicates global nodal forces and f element nodal forces. Making use of Equation 


B.1-4, the equilibrium equations Equation B.1-5 follows, 
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elel 





F=f = ku, —ku,, 
elel ele2 

F=f +f, = —ku,+(k,+k,)u, -—k,u,, B.1-6 
ele2 

F =f; =—k,u, +k,u, 


In matrix form 





B.1-7 


The two individual stiffness matrices for element 1 and element 2 are indicated in Equation 


B.1-7. Equation B.1-7 is the global set of equations, with K the global stiffness matrix, D 


the global nodal displacements and F the global nodal force vector. 


An alternative way of assembling the global stiffness matrix follows: The stiffness matrix for 


each element can be “enlarged”, 


-k k Ow |=1F, for element 1 


| 0 0 Olu 0 
0 0 0 | u, 0 
0 k —-k,|u, |= oe for element 2 B.1-8 





The two equations can now simply be added together, 
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k, -k 0 u fi F, 

-k ktk, -k, |u, |=| f +f |=|F, B.1-9 
ele2 

0 -k, k, || us Í F, 


B.1.3 Boundary Conditions 


Can Equation B.1-9 be solved for the nodal displacements? In order to accomplish this, we 


need to calculate the inverse of the stiffness matrix, D = K R 
The determinant of the stiffness matrix is given by, 
detK = k, (k +k, k, -k |- kjk, = kk, +k,-k,-k,k,=0 B.1-10 


The determinant is equal to zero, which means that the inverse of K is not defined and the 


nodal displacements can not be solved for. 


In physical terms this means that the system (springs) are not constrained, i.e. rigid body 
motion is possible. To avoid rigid body motion, constraints should be applied. The most 


common constraint is to “fix” one or more of the nodes. 


For this example, let’s fix node 1: u; = 0 and F = F; = P. The external force applied at 
node | would be the reaction force, an unknown. Equation B.1-7 can be written as 

k —k, 0 [o0] TF 

-k k+k, -k, |u, |=| P 

| 0 -k, k, || u, P 


which reduces to 


[ki +k, =k, |u, P 
= B.1-11 
[Lees k, || u, P 
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and F, =—k,u,, the reaction force. This is equivalent to deleting row 1 and column 1 in 


the stiffness matrix. If it was node 2 that was constrained, it would have been row 2 and 


column 2 that would have been deleted. 


The reduced stiffness matrix is no longer singular, and the system can be solved for the 


unknown nodal displacements at node 2 and node 3 (node | is fixed). 
1 1 2P 
~ | 1/ 41 | 2P/ 4P ne 
Te, Jat Ne Jii 7 


and the reaction force follows from Equation B.1-11, F, =—k,u,=-2P, which ensures 


u3 u3 


equilibrium. 


Note: The unconstrained system has 3 degrees-of-freedom (one displacement per node). The 


constrained system has only 2 degrees-of-freedom. 


Now, with the knowledge of assembling and boundary conditions know, we can move onto 
more complex elements. The process of assembling a global stiffness matrix and applying 
boundary conditions are the same for any type of element, it is just the formulation (and size) 


of the stiffness matrix that may differ. 


2.2 Stiffness Matrix Formulation: Bar Element 


The bar element is similar to a spring element, i.e. it can only resist an axial force. Each bar 
element has 2 nodes, i.e. 2 degrees-of-freedom. Let’s look at a simple bar element with 


constant area A, Young’s modulus E and length L: 


o F 
o=Eée > €é=-—, o=— > §€=— 
E A AE 
but also 
AL FL AE 


é=— > AL=el > AL=— > F=—AL 
L AE L 


Comparing this to a spring element, 
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F = KAL + k= B.2-1 


2.2.1 The Direct Method 


The direct method of calculating the stiffness matrix was used for the spring element. Above 
it has been shown that the stiffness of a bar element is given by Equation B.2-1. Making use 


of Equation B.1-2, the stiffness for the bar element follows as, 


1 -1 
— aa 2.2-2 
L-1 1 


This, of course, assumes small deformations and elastic materials. 


2.2.1 The Formal Procedure 





u = Nju] + Noup 








Figure 2.2-2 — Shape functions N, and N, for a two-node 


Let’s write an expression for the axial displacement u of an arbitrary point on the bar. 

Assume that node 2 is fixed, and node 1 is given a unit displacement, Figure 2.2-2a. The 

displacement at node 2 will be zero, and the displacement at node 1 should be unity (1). The 

deformation, for example, in the centre of the bar should be 0.5. This is written as 

_L-x 
L 


u 





Mi B.2-2 


where x starts from node | as in the figure, and u, is the displacement of node 1. The same 


can be done with node 1 fixed and node 2 given a unit displacement, 
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u= fy, B.2-3 
L 


Combining Equation B.2-2 and B.2-3 in matrix notation 


= u E 
u= pee oli or u=N d 2.2-5 
E Lll = 


where N is called the shape function matrix and d the vector of element nodal degrees-of- 





freedom (translation in this case, but could also be rotations as in beam elements). 


Let’s now return to the definition of strain given in Chapter 1 (notes). Since this is a one- 


dimensional case (only axial deformation), there is only one strain component namely 





gpa B.2-4 
Ox 
With u given by Equation 2.2-5, 
2-2 a] 
ôx Ox 
but d is the nodal displacements (constants), and not a function of x, 
pe A eng 2.2-6 
ôx | Ox 


where the matrix B is called the strain-displacement matrix, containing shape function 


: 5 1 ; 
gradients, B = |- T a . The bar element has constant strain. 


To formulate the stiffness matrix, we need to look at energy principles (see section 3.1 of the 
textbook). First assume a linear elastic spring/bar, with the force-displacement behavior as in 


Figure B.2-1. 
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Figure B.2-1— Force —displacement relation for simple spring/bar element 


The increment in work done on the bar is given by 


dW = F du 
but F=ku + dW =kudu >W = fkudu= Shu =F B.2-5 


This is the total work done on the bar when extended by a distance u. In the case of multi 


degrees-of-freedom, the work is given by 

F B.2-6 
where d is the vector containing nodal displacements u. The work done on the bar, Equation 
B.2-6, goes into increasing the strain energy. The increment in strain energy per unit volume 


is given by 

dE, =ode 
but o =E¢=EBd (using Equation 2.2-6) 
thus £, =| odé= [ze dé = Le 

2 

In the case of multi degrees-of-freedom, this becomes 

Bo=+5" Re 

2 


The total strain energy is calculated over the total volume of the element, 
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E=/ Leea B.2-7 
12 


Making use of Equation 2.2-6 


i B.2-8 
zeg | B'EB wr] d 
2 y 


Conservation of energy states that the work (Equation B.2-6) input should equal the change in 


strain energy (Equation B.2-8), 





W =E 
ae B.2-9 
sag Beate [BEB dy | d 
2 2 7 
From this it is clear that 
F= ( BEB wr) d B.2-10 
V 
which is similar to the one-dimensional expression F = k d, and we conclude that 
k=|B'EB av 2.2-4 
V 


which is the general formulation for the stiffness matrix of an element. 


Let’s now apply this formulation to calculate the stiffness of the bar element. Using the 


definition of B from Equation 2.2-6, 


T 
E-{|-4 L| £|-4 | aa 
sL L L L L 


where the cross section area A is taken to be constant. 
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1 toot toot 
a tle- L Ads = E Al L L dx= AEL L L 
L É E É E 
2.2-7 
a. 
-apl L i -4 7 
eee eee a Ee A 
L L 


This result is the same as obtained using the direct method, Equation 2.2-2. 


Note: The stiffness matrix formulation given in Equation 2.2-4 can be applied to any element. 
For different elements it would only be the strain-displacement matrix B that differs and for 


two-dimensional and three-dimensional problems Æ is replaced by a matrix containing both £ 


and Poisson’s ratio v. 
The process of assembly and constraints (reduced stiffness matrix) developed for spring 


elements can also be applied to the bar elements developed in this section. It is just the 


formulation of the element stiffness matrix that is different (although the same size 2 x 2). 
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2.3 Stiffness Matrix Formulation: Beam Element 


2.3.1 The Direct Method: Simple Plane Beam 


Figure 2.3-la shows a simple plane beam element. The elastic modulus Æ and moment of 
inertia J are constant along the length of the beam. Only lateral displacements are assumed v 
= v(x). For a beam only loaded at the end nodes with forces and moments, it can be shown 
from elementary beam theory that the lateral displacement is a cubic function in x. Each of 


the two nodes has three degrees-of-freedom, two translational and one rotational. 

















Figure 2.3-1 — Simple plane beam element and its nodal dof and shape function definitions 
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When the stiffness matrix for the bar element was developed, it was shown that the j™ column 
in the stiffness matrix represents the nodal forces (at each of the nodes respectively) that 
should be applied to have unit displacement (rotation) at the  degree-of-freedom and zero 
displacement (rotation) at all the other degrees-of-freedom. This procedure will be followed 


to calculate the beam stiffness matrix. Assume the degrees-of-freedom are assembled as 


d=|y, A Vv, aie 


The first degree-of-freedom is vı. Set vı = 1 and all other degrees-of-freedom equal to zero. 
This is indicated in Figure 2.3-1c. Forces and moments should be applied at all degrees-of- 
freedom. What should these forces and moment be? Let’s return to Basic beam theory” as 


indicated in Figure B.3-1. 


Slope Deflection 
m j MI? 
A a we MI at B eT ae B 
LV I 2EI 
——_— l 'B 











Figure B.3-1— Basic beam theory 


Using Figure B.3-1, the lateral displacement of node | (Figure 2.3-1c) can be written in terms 


of the applied moment and force at node 1: 





3 2 
pl fat B.3-1 
3EI 2EI 
the slope (angle) at node 1 should also be zero: 
2 
6, je B.3-2 
2EI EI 





> Benham, PP, Crawford, RJ, Armstrong, CG, Mechanics of Engineering Materials, Second Edition, 1996, 
Addison Wesley Longman Limited. 
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Solving these two equations simultaneously for kı; and ky; 


12ET 
ki = 3 
L 





6ET 
kz = 





2 


L 


From equilibrium, we can write the following two equations, 


Sum of forces in y-direction: 


Sum of moments around node 2: 





Using Equation B.3-3, these can be solved: 


_12E1 
i 





ky = 





0= ki T kz 
0= kz S ky z kL 
6ET 
LC 


B.3-3 


B.3-4 


B.3-5 


Following this procedure for all the other degrees-of-freedom, the final element stiffness 


matrix is obtained, 














[ 121 6EI 12EI 6EI 
T r j v 
6E1 4E1 6EI 2E1 
v L Lv L 
Koan = EL 6EI NEI 6EI 
LC Ï LC i 
6E1 251 6EI 4E1 
JA L Te L 





Note: This element has 4 degrees-of-freedom, two per node. 


The equilibrium equation for this element 

















EI 6E1 EI 6EI 
if Lv Lv r 
6EI 4E1 6EI 2E1 
Te L T L 
EI 6EI REI 6EI 
T T T r 
6EI 2E1 6E1I 4E1 




















B.3-6 


15 of 27 


Eindige Element Metodes 414 / Finite Element Methods 414 


2.3.2 The Formal Procedure: Simple Plane Beam 


In the case of a beam element the stiffness matrix is given by an expression similar to that in 


Equation 2.2-4, 


k= f B EIB dV 2.3-3 
V 


where Bis now the matrix given the relation between nodal displacement and beam 


ð f a at : : 
curvature A The lateral displacement of the beam is given in terms of generalized 
ox 


coordinates 2, 


v= B+ Byxt+ Bx + B,x 2.3-4 


The f’s can be written in terms of nodal displacements/rotations making use of a similar 








procedure used for bar elements. For example at x = 0, v= vı and 0=60, = z ; 
x x=0 
The result is given by 
a 
0, z 
v=|[N, N, N, N,] =Nd 2.3-5 
Vo 
L8] 
The shape functions are given in Figure 2.3-1. The curvature of the beam element, 
2 2 
OV! gee 2.3-6 
dx dx 


where the strain-displacement matrix is (using the shape function definitions) 


2-| 6 12x 4 6x 6 12x nS 








2.3-7 


E A a ee 
After substitution of Equation 2.3-7 into 2.3-3, and some manipulation, the same stiffness 


matrix as in Equation 2.3-2 is obtained. 
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This formulation gives an exact solution to the beam problem under the following conditions: 
e The beam is initially straight 
e Linear elasticity assumed 
e No taper — constant cross section 


e Under small deformations 


; er M: ; 
The stress in the beam is given by the well-known formula o = T where the moment M is 


given by 
of 
M =EI—~=EIBd 2.3-8 
dx 


2.3.2 2D Beam 


A 2D beam is a combination of a bar element and a simple plane beam element. It can resist 
axial loads, transverse loads and bending moments. The equilibrium equation (and stiffness 
matrix) is a combination of that of the bar element (Equation 2.2-7) and that of the simple 


plane beam (Equation B.3-6) 

















[ AE _ AE ilies) VP 
L 0 0 L 0 0 u Pran 
12EI 6EI _ 12EI 6EL 
0 Cc E 0 ig E vı Fani 
6EI 4EI _ GEL 2EI 
0 FE TE 0 P 7 0, M, 
ini pi = B.3-7 
-7 0 0 + 0 O lft.) | Faia 
_ 12EI _ 6EI 12EI 6EI 
0 Fo E 0 E lA v 2 trans2 
6EI 2EI _ 6EI 4EI 
n To o ro T |a] LM, | 
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2.4 Properties of the Stiffness Matrix: Avoiding Singularity 


The element stiffness matrix k and the global stiffness matrix K are symmetric. This is, 


however, only true if the material behaviour is linear. The diagonal components of both 


stiffness matrices are always positive. 


If a structure is unsupported or inadequately supported, K can be singular, and the system 
equations unable to be solved. To prevent this, the structure must be supported (constrained) 
to avoid rigid body motion. A structure may also have a singular K if it contains a 


mechanism. 
2.5 Mechanical Loads: Stresses 


Loads can be either applied forces or moments at a point, or surface pressure, or body forces 
(inertia, gravity). Concentrated forces and moment are applied at nodes (point loads). 
Moments, however, can only be applied at the node if at least one element connected to this 


node, has a rotational degree of freedom, e.g. a beam element. 


Distributed loads, such as pressure, acts on elements between the nodes. These distributed 
loads must be converted to equivalent nodal loads, since in FEM, loads can only be applied at 
nodes. For bar elements, the following procedure can be used to obtain work-equivalent 


loads. 


Take a distributed load q [N/m], acting on a bar element as indicated in Figure 2.5-la. The 
load q causes the bar to extend. The nodal displacements would be u; and u2. The work done 
by the distributed load is given by force times displacement. With the nodal displacements 


known, the displacements within the element follow from the shape functions, u(x)=d N. 


The work is given by fugo) dx = fa N q(x) dx. We are now looking for nodal forces F, 
L L 


that would result in the same amount of work being done - work-equivalent nodal loads. 
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d F, = [aN qe) dx B.5-1 
L 


but since d are constant nodal displacements 


aF, =a[| Nao a| > F,=|N q(x) dx B.5-2 


In the case of a point load F, it becomes 


d F,=|dNF dk >F,=NF B.5-3 
L 


Element Structure 


i q 2 
pee de: pee eer 


een ee ne be .,, —_______>| 





Figure 2.5-1 — Uniformed distributed axial force q on a two node bar element. 


The following example demonstrates the procedure where a distributed and a point load is 
applied to a single bar element, Figure B.5-1. The distributed load q varies linearly from q at 


node 1 to qo to at node 2, 


L-x x 
xy =*= q +> B.5-4 
q(x) TS 


l ; . : 2L 
and a point load F is applied at a distance x = a The work equivalent nodal forces are 


calculated using Equation B.5-2 and B.5-3 


F,=(N q(x) dx +N F B.5-5 
L 
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q2 








Figure B.5-1 — Equivalent nodal loads for a single bar element 


L-x 





Making use of the shape function definition of a bar element, Equation 2.2-5 N = i ; 
L 


Equation B.5-5 becomes 




















[L-x] L-x 
z. L L-x xX i L 
F=] x í z qı Eq, Ja x F 
LL] L L 2L 
3 
[L-x] [L-4 
_ L 7% x L 
F=] x 7 nta a ab F 
L | LOOL 











2q, +q F 
| l /s B.5-6 


F, == + 
qı +2q, a 


6 
In the case where the load is uniformly distributed, q; = q2 = q, and F = 0, 


wp 
F=- £2 B.5-7 


This is demonstrated in Figure 2.5-1b, where the total force is gl. Figure 2.5-1c and d shows 


how the nodal loads add together when more than one element is assembled. 
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2.6 Thermal Loads: Stresses 


When the temperature of an unrestrained body is uniformly increased, the body will deform, 
but the stresses will be unchanged. Temperature gradients are, however, more complicated, 
and thermal stresses can occur even if the body is unrestrained. In this course only a uniform 
temperature change will be treated. Software does these calculations automatically, with the 
user only specifying the temperature change. The next example, however, shows how the 


thermal changes are handled. 





Figure 2.5-1 — Two element bar model loaded by externally applied load P and by uniform heating an 


amount AT 


Figure 2.5-la shows two bar elements assembled together. A load P is applied at node 2 in 
the negative x-direction, and a load P at node 3 in the positive x-direction. Let the cross 
sectional area A, elasticity Æ and the coefficient of thermal expansion a be constant over the 


bar. 


Due to the temperature change alone, the bar will extend, but the stress will be unchanged. 
The change in strain can be modelled by applying nodal loads, but these loads will also cause 
the stress to increase. The solution is to apply thermal loads, that should yield the same strain 
deformation as the temperature change, and then to subtract the stress caused by the thermal 
load. 


The definition of the coefficient of thermal extension ¢=a@AT is used to calculate the 


thermal loads, 


F, =0A=EcA=a AEAT B.6-1 
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This load would cause the same increase/decrease in strain as the increase/decrease in 


temperature AT . The element stiffness matrix is given by 


a. _ Ly at ~! es - F; 
k, =k, =— and the thermal load r,, = 7, = 2.6-1 
L\-1 1 F, 
Assembly of elements yields the global equilibrium equation 
1 -1 0 || —F, +R, 
£ -1 2 He apa HP 2.6-2 
0 -1 1 jiu, Er FP 


where R, is the support reaction at node 1. The stiffness matrix in Equation 2.6-2 is singular. 


Singularity is removed by the support conditions, i.e. removing row 1 and column 1. 


2 —=—lļ|u -P 
AE Pies 2.6-3 
L\-1 1 |lu,| |F,+P 


Solving this for the nodal displacements, and making use of Equation B.6-1, 


PL 
u,=0, u,=a@LAT, u, =2aLAT +— 2.6-4 
AE 
To calculate the stress in each of the elements, we have to remove the stress due to thermal 
loads. 
U, =U 
on = Ee, -aEAT =E -aEAT =0 2.6-6a 
uU, — U, P 
Oy, = He, — aE AT =E 2 aE AT == 2.6-6b 


These results can be checked through inspection! 
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2.7 Transformations 


The textbook only handles transformations in Chapter 4, section 4.3. However, it fits better 


in Chapter 2 and is needed before a complete bar/beam example can be done. 


We have to diffirentiate between local coordinates and global coordinates. The user defines 
the geometry of a FEM model in the global coordinate system XYZ. Software typically 
generates an element stiffness matrix in local coordinates xyz, then automatically converts to 
global coordinates for assembly of elements. Global and local systems may be parallel or 


even coincident. 


The stiffness matrix of a element is most easily written in local coordinates, but it may be 
arbitrary orientated in global coordinates. Rather than formulate the element stiffness matrix 
in global coordinates, it is easier to transform the stiffness matrix to global coordinates. 
Transformation is carried out outomatically by software. Let’s first look at the transformation 


of a single point in two-dimensional space. 








Figure B.7-1— Coordinate transformation of a single point in two-dimensional space 
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From figure B.7-1, the following relations can be written, 


x, =X, cos0+Y, sind 


; B.7-1 
Y, =Y, cosO—X, sin 
In matrix notation 
xX, cos sin@ || X, 
|| ly B.7-2 


ie oe 


Pp 


where T is called the tranformation matrix which is orthogonal: T 


2.7.1 Bar Element 


Let u’; and u’? be the nodal displacements in local coordinates of a bar element with two 
nodes and (u;, v;) and (u2, v2) the nodal displacements in global coordinates. The 


transformation becomes 








u 
ui c s 0 Oliv 
= B.7-3 
u, 0 0 c s||u, 
L¥2 
or d'=Td where c=cos0,s=sin@ and T is build up from components in T ; 


Equation B.7-2. Similarly, the nodal forces can be written as 


S'=Tf B.7-4 
The equilibrium equation in local coordinates 

AE| | oiu) | fa 

Piles ADI d 

B.7-5 

k d'=f' 
Using the transformations defined in Equation B.7-3 and B.7-4 

k Td=Tf B.7-6 
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B.7-7 


where the force f is in global coordinates and the stiffness matrix in global coordinates is 


given by 
is 2 
c cs 
Ee cs s 
Teren 
Bagn mes 
2 
-ces -s 





2 
—Cc —cs 
—-cs —S 
2 
Cc cs 
2 
cs S 


The element stress is calculated using 


o- Be =EBd= 5) 

















B.7-8 


B.7-9 
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2.7.1 2D Beam Element 


The equilibrium equation for a 2D beam element is given in Equation B.3-7, 























T 0 0 =T 0 0 [u | | Fian 
12EI 6EI IAEI GEI ' ' 

0 P 2 0 T 2 vi transl 
6EI AEI _6EL 2E ' 

0 %4 4 2 jjo M! 

= B.7-10 
A 0 0 T 0 0 us FB ois 
12EI 6EI 12EI 6EI 

0 E P Te 0 P ~ P v, kana 
SEI 2EI — SEIL 4EI ' ' 

0 g u 9 et a. | Mi | 


Both the u and v displacements need to be rotated in the case of a local coordinate system not 
aligned with the global coordinate system. The rotational degree-of-freedom ©, however, 
remains unchanged since the local and global z-axis (out-of plane) remains aligned no matter 


what the rotation (2D rotation). For a beam element, the transformation matrix can be written 


as 

lc s 0 0 00] 
-s c 0 0 0 0 
_ 0 01 0 0 0 

T= B.7-11 
0 00 c s0 
0 0 0 -s c0 
0 00 0 01 


which is build up from components in ia , Equation B.7-2. 
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Eindige Element Metodes 414 / Finite Element Methods 414 


The following example demonstrates the use of coordinate transformation for bar elements. 


Element 1 





Figure B.7-2 — Example: Coordinate transformation 


Figure B.7-2 shows two identical bar elements: A, E, length = L. 


Find: 


1) The unrotated stiffness matrix of each element. 
2) The rotated stiffness matrix of each element. 

3) The global unconstrained stiffness matrix. 

4) The global constrained (reduced) stiffness matrix. 
5) The global force vector. 

6) The displacement of node 2. 


7) The stress in each bar. 
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e CHAPTER 2 Finite Element Methods 414 
Example: Coordinate Transformation 


kPa := 1-10" Pa kN := 1-10°N 


MPa := 1.10°Pa MN := 1-10°N 


PROBLEM: 


Pı 





Element 1 








Find 1) The unrotated stiffness matrix of each element. 
2) The rotated stiffness matrix of each element 
3) The global unconstrained stiffness matrix 
4) The global constrained (reduced) stiffness matrix 
5) The global force vector 
6) The displacement of node 2 


7) The stress in each bar. 
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Solution 1: 


In local coordinates the stiffness matrix is: 


Uy: ND 
peep _AE 1 -1) uy 
ica aa a ie 


Solution 2: 


For element | the rotation angle is 6, := 45deg 


cos(0} = 0.707 2 = 0.707 
2 


Using Equation B.7-8 (Lecture notes) 


uj V1 "2 Y2 Global numbers: 


For element 2 the rotation angle is 0, := 135deg 


uz, S28 98 

1-11-11) w 
AE|-!1 1 1 -1 Y2 
2L|-1 1 1 -1 

i aii A vy 


Solution 3: 


1 1 -1-1 0 0% F1 
1 1 -1-1 0 0|] Fy 
AE|-1 -1 2 0 -1 1 ||% Fy 
pate E es ed ce 
0 0-11 1 -1{ly] | Fe 
0 0 1 -1-1 1Jiv3) (F 
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Solution 4: 


Node 1 and Node 3 are fully constrained in the x- and y-directions: Delete row and column 
and 5,6 respectively: 


Moy 
K, AE k; 0) % 
duced = 77 
educe: 2L l0 2) 2 
Solution 5: 
P; ) 
Freduced = Py J 


Solution 6: 


e é ar 
vy) [2L (0 2) 


o i 
v2) ~ AE Po) 


Solution 7: 


fa 
Py J 


Using Equation B.7-9 (Lecture notes) 


0 ) 


(Py + P2) for element 1 





(Py = P2) for element 2 
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